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Abstract 

> 

In 2009, for the first time in Cape Verde, an outbreak of dengue was reported and over 
twenty thousand people were infected. Only a few prophylactic measures were taken. The 
effects of vector control on disease spreading, such as insecticide (larvicide and adulticide) 
and mechanical control, as well as an hypothetical vaccine, are estimated through simulations 
with the Cape Verde data. 

Keywords: dengue; Cape Verde; vaccine; insecticide; mechanical control; basic reproduction 
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j_] 1 Introduction 



Dengue is a disease which is now endemic in Africa, America, Asia, and the Western Pacific. In 
Europe there has been no registered cases so far, yet the main vector of the disease is already 
present there and has been followed in Madeira. According to the World Health Organization 
www. who . int/topics/dengue/en| 2011), the incidence of dengue has grown drastically in the 



last decades and roughly two fifths of the world population is now at risk. So far, the strategies 
to control Aedes aegypti mosquito proved to be inefficient and the continuous level of surveillance 
of the mosquito is low. Climatic changes, unsanitary habitat, poverty, uncontrolled urbanization, 
and global travel favor the propagation of dengue fever. Inadequate water supply requires large- 
scale water storages, which are ideal breeding habitats for mosquitoes. The increasing movement 
of people and goods has enabled the dengue virus and its vectors to spread to new parts of the 



world (Semenza and Menne 2009). There is no effective vector control and dengue infection rates 



have been increasing during the last 40 years, e.g., in Thailand has increased from 9/100,000 in 
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1958 to 189/100,000 in 1998 ( |Healstead[ |1992[ |Endy et al.[|2002] . The virus emerged for the first 
time in Cape Verde at the end of September 2009. The outbreak was the biggest ever recorded 
in Africa. As the population had never had contact with the virus, herd immunity was very low. 
Dengue type 3 spread throughout four of the nine islands. The worst outbreak occurred on the 
island of Santiago, where most people live, up to 1000 cases per day in November. The Ministry 
of Health of Cape Verde reported over 20,000 cases of dengue fever, which is about 5% of the 
total population of the country, between October and December 2009. From 173 cases of dengue 
haemorrhagic fever reported, six people died ( www . dengue . gov . cv( 2011 [www . minsaude . gov . cv[ 
20TT|. 



2 Dengue and Aedes aegypti 



Dengue is transmitted to humans by mosquitoes, mainly Aedes aegypti, and exists either as dengue 
fever (DF) or as dengue hemorrhagic fever (DHF). The disease can be contracted by one of four 
types of viruses. Infection with one serotype confers lifetime immunity against that serotype, but 
not to the others, and there is evidence that a prior infection increases the risk of developing 
DHF for people infected with another serotype. DF is characterized by sudden high fever (3 to 
7 days) without respiratory symptoms, intense headache, and painful joints and muscles. The 
haemorrhagic form is also characterized by sudden fever, nausea, vomiting, and fainting due to 
low blood pressure by fluid leakage. It usually lasts between two and three days and can lead to 
death. There is no specific effective treatment for dengue. Fluid replacement therapy is used if 
clinical diagnosis is made early. Vaccine candidates are undergoing clinical trials ( www .who . int/^ 



topics/dengue/en 2011). 

Aedes aegypti is closely associated with humans and their dwellings, thriving in crowded cities 
and biting primarily during day light. Humans also provide nutrients needed for mosquitoes to 
reproduce through water-holding containers, in and around human homes. In urban areas, Aedes 
mosquitoes breed on water collections in artificial containers such as plastic cups, used tires, 
broken bottles, and flower pots. With urbanization, crowded cities, poor sanitation, and lack of 
hygiene, environmental conditions foster the spread of the disease which, even in the absence of 
fatal forms, breeds significant economic and social costs (absenteeism, immobilization, debilitation, 



and medication) ( Derouich and Boutayeb 2006 Ferreira et al.[ 2006 ) 



Female mosquitoes acquire infection by taking a blood meal from an infected human. These 
infected mosquitoes pass the disease to susceptible humans. Female mosquitoes lay their eggs on 
inner wet walls of containers. Larvae hatch when water inundates eggs. In the following days, 
the larvae feed on microorganisms and particulate organic matter. When the larva has acquired 
enough energy and size, metamorphosis changes the larva into a pupa. Pupae do not feed: they 
just change in form until the adult body, the fling mosquito, is formed. The newly formed adult 
emerges from the water after breaking the pupal skin. The entire life cycle, from the aquatic phase 
(eggs, larvae, pupae) to the adult phase, lasts from 8 to 10 days at room temperature, depending 



on the level of feeding ( Christophers , 2009 1 



It is very difficult to control or eliminate Aedes aegypti mosquitoes because they quickly adapt 
to changes in the environment and they quickly pullulate again after droughts or prophylactic 
measures. The transmission thresholds and the extent of dengue transmission are determined by 
the level of herd immunity in the human population to circulating virus serotype(s), virulence char- 
acteristics of the viral strain, survival, feeding behavior, abundance of Aedes aegypti, climate and 
human density, distribution, and movement (Scott and Morrison 2004). There are two primary 



preventions: larval control and adult mosquito control, depending on the intended target (Natal 



2002). Larvicide treatment is an effective control of the vector larvae, together with mechanical 



control, which is related to educational campaigns to remove still water from domestic recipients 
and eliminating possible breeding sites. The larvicide should be long-lasting and preferably have 
World Health Organization clearance for use in drinking water (Derouich et al. 2003). The ap- 



plication of adulticides can drop the mosquito vector population. However, the efficacy is often 



constrained by the difficulty in achieving sufficiently high coverage of resting surfaces (Devine 
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et al. 2009). 



While [Feng and Velasco-Hernandez ( 1997 ) investigate the competitive exclusion principle in a 
two-strain dengue model, Chowell et al. ([2007 ) estimate the basic reproduction number for dengue 
using spatial epidemic data. Tewa et al. (20091, established the global asymptotic stability of the 
equilibria of a single-strain dengue model. Thome et al. (2010) introduce a sterile insect technique. 



3 Model 



We adapt the model presented in Dumont et al. (2008) and Dumont and Chiroleu (2010) to 
dengue, with the considerations of Rodrigues et al. (2009 2010a|b |. We consider three controls 
simultaneously: larvicide, adulticide, and mechanical control, with mutually-exclusive compart- 
ments, to study the outbreak of 2009 in Cape Verde and improve upon IRodrigues et al. (2009). 
We denote Sh the total number of susceptible, //j the total number of infected and infectious, 
and Rh the total number of resistant individuals. The total human population is a constant 
Nh = Sh{t) + Ihit) + Rh{t)- The population is homogeneous, which means that every individual of 
a compartment is homogeneously mixed with the other individuals. Immigration and emigration 
are ignored. 

Female mosquitoes are in total number in the aquatic phase (including egg, larva, and 
pupa stages), Sm are susceptible, and are infected. Mosquitoes are considered to live too 
briefly to develop resistance. Each mosquito has an equal probability to bite any host. Humans 
and mosquitoes are assumed to be born susceptible. 

We consider three controls: the proportion ca of larvicide, the proportion Cm of adulticide, and 
the proportion a of mechanical control. Larval control targets the immature mosquitoes living 
in water before they bite. The natural soil bacterium Bacillus thuringiensis israelensis (Bti) is 
sprayed from the ground or by air to larval habitats. The control of adult mosquitoes is necessary 
when mosquito populations cannot be treated in their larval stage. Depending upon the size of 
the area, either trucks for ground adulticide treatments or aircraft for aerial adulticide treatments 
are used. The purpose of mechanical control is to reduce larval habitat areas. The parameters 
used in the model are: 

Nfi total population; 

B average daily biting (per day); 

Prnh transmission probability from (per bite); 

Phm transmission probability from (per bite); 

l//ift, average lifespan of humans (in days); 

l/i]fi mean viremic period (in days); 

l//i„i average lifespan of adult mosquitoes (in days); 

(p number of eggs at each deposit per capita (per day) ; 

1/ fJ-A natural mortality of larvae (per day); 

rjA maturation rate from larvae to adult (per day); 

m female mosquitoes per human; 

k total number of larvae per human. 

The dengue epidemic is modelled by the nonlinear time-varying state equations 



S'f^it) - ^ihN,, - + ^ih) Shit) 



R'^it) = rihlh{t) - lihRh{t) 



(1) 



and 



S'mit) 



BPhm pf ' Sm {t) — (/irn 



ImXt)) - iVA + (J-A+Ca) Am(t) 

Sm{t) 



(2) 
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Figure 1: Epidemiological model SIR (Susceptible, Infected, Recovered) + ASI (Aquatic phase. 
Susceptible, Infected) with Sh, Ih, Rh, ^m, Sm, and /,„ as state variables, ca, Cm, and a as 
controls, and Nh, B, P^h-, Phm, fJ-h, Vh, fJ'm, ^, fJ-A, Va, m, and k as given parameters. 



with the initial conditions 



S'h(O) — 5/10, lh{0) — IfiO, Rh{0) — RhO, 
^m(O) = A,nQ, S„i{0) — SmO, ^m(O) = ImO- 



Figure [T] shows a scheme of the model. 



4 Equilibrium points and the basic reproduction number 

We study the solutions of System ([T])-([2]) in the closed set 

n = {{ShJh, Rh,Ara, Sm, Im) G ^\ ■ Sh + h + Rh < Nh, Am < kN,,, Sm + I.n < mNh} . 

The n set is positively invariant with respect to Eq. ([T])-([2| ( |Rodrigues et al. 2012). System 
([l])-([2]) has at most three biologically meaningful equilibrium points (cf. Theorem [2|) . 

Definition 1. A sextuple E — {Sh, Ih, Rh, Am, Sm, Im) G is an equilibrium point for System 
if it satisfies 

fJ-hNh - (^Bl3mhj^ + fJ-h^ Sh^O 
BjSmhj^Sh - [rjh + ^J.h)Ih = 

^hIh - iJ'hRh = 

^{}~ 5¥f;:) + Im) - {flA + A*A + CA)Ar, 

r]AAm - i^BPhm + Mm + 



= 



(3) 



s„ 



An equilibrium point E is biologically meaningful if and only if _E G O. The biologically meaningful 
equilibrium points are said to be disease-free or endemic depending on //, and /,„ : if there is no 
disease for both populations of humans and mosquitoes {Ih ~ Im = 0), then the equilibrium point 
is a disease-free equilibrium (DFE); otherwise, if //j > or /,„ > 0, then the equilibrium point is 
called endemic. 

Tlieorem 2. System ([l])-([2]) admits at most three biologically meaningful equilibrium points, at 
most two DFE points, and at most one endemic equilibrium point. Let 



fHA) . 



(4) 
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and Et 



El = (7V^, 0,0, 0,0,0) 

m rn 



Af/,.,0,0 



X = akB PhniPnihM, 
akNhM akNhM 



VA^ ' (Mm + 



,0 



(5) 
(6) 



S 



It. 



BPhmi-akBPmhM - 



n 

+ c„0) 



+ Cm) 



+ ^ih)Bl3h„ii-akBl3„^hM - (pHhijJ-r, 
VhNh - x) 



NhkaM 



{-akBl3mhM - ipiXhilJ-r, 



.))' 



0)' 



(7) 



iV^afc {cmiph + 'r)h)ipA + 'riA+ ca) + A^n 



m rjh + HhBPh 
A^fta/c {-r]ip{^j.h + Vh) + A^m(??/iMvi + M^ca)) 



{^pBPhm^J■h{^J■m + Crn) + ' 

^ 0; then there is only one biologically meaningful equilibrium point E^, which is a DFE 
point. If A4 > with S, > Xj then there are two biologically meaningful equilibrium points E^ 
and E2, both DFE points. If Ai > with ^ < X> then there are three biologically meaningful 
equilibrium points E^, E2, and E^, where E^ and E2 are DFEs and E^ is endemic. 

Proof. System ([s]) has four solutions obtained with the software Maple: E^, E2, E^, and 
The equilibrium point E^ is always a DFE because it always belong to O with 1}^ = 1^ = 0. In 
contrast. El is never biologically realistic because it always has some negative coordinates. The 
other two equilibrium points, E2 and E^, are biologically realistic only for certain values of the 
parameters. The equilibrium E2 is biologically realistic if and only if > 0, in which case it is 
a DFE. For the condition < 0, the third equilibrium E^ is not biologically realistic. If > 0, 
then three situations can occur with respect to £"3: if ^ = x? then E^ degenerates into E2, which 
means that E'^ is the DFE £'2; if ^ > Xi then E^ is not biologically realistic; otherwise, one has 

□ 



£3 G Q with //i 7^ and /,„ 7^ 0, which means that E^ is an endemic equilibrium point. 



By algebraic manipulation, AI > is equivalent to the condition 



{flA + MA + CA){p-rn + Cm) 



> 



fVA 



1, which is related to the number for offspring mosquitos. Thus, if Ai < 0, then the mosquito 
population will collapse and the only equilibrium for the whole system is the trivial DFE E'^. If 
A^ > 0, then the mosquito population is sustainable. From a biological standpoint, the equilibrium 
E2 is more plausible, because the mosquito is in its habitat, but without the disease. 

Definition 3. (Hethcote, 2008) The basic reproduction number, denoted by TZq, is defined as 
the average number of secondary infections occurring when one infective is introduced into a 
completely susceptible population. 

The basic reproduction number provides an invasion criterion for the initial spread of the virus 
in a susceptible population. For this case. 



Theorem 4. The basic reproduction number TZq associated to the differential System ([I]) -([2]) 



IS 



7^n = 



akB'^ PhmPmhM 



(8) 
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Proof. In agreement with van den Driessche and Watmough ( 2002 ) , we consider the epidemiolog- 



ical compartments with new infections, Ih and Im- The two differential equations related with 
these two compartments are rewritten as x'(t) = !F{x{t)) — V{x{t)), where x = {Ih, Im), ^ is the 
rate of production of new infections, and V the transition rates between states: 



The Jacobian derivatives are 

Bp„,h^ \ _ f Vh+^^h 



T-1 



The quantity Jj^{x}Jy(^x^ gives the total number of new infections over the course of an outbreak. 
The largest eigenvalue gives the asymptotic growth of the infected population, giving TZq as the 
spectral radius of the matrix Jj^{x)J^^x) ^ DFE point. Maple was used to obtain 



The basic reproduction number TZq in Eq. ^ is obtained from replacing Sh^pE ■^'todfe 
Eq. (|9| by those of the DFE E^. □ 

The model plays on the populations of host and vector, and the expected basic reproduction 
number should reflect the infection transmitted from host to vector and vice-versa. Accordingly, 
TZo can be seen as TZq = {TZhm x TZ„ih)^ ■ The infection host-vector is represented by TZhm = 
fiffhrng mpFE ^ -v^rherc the term •^'^'""^'"dfe represents the transmission probability of the disease 



from humans to mosquitoes in a susceptible population of vectors, and the term ^^^^^ the human 
viremic period. Analogously, the infection vector-host is represented by TZmh — ^^^'(c^+i^'^^) ' where 



represents the transmission of the disease from mosquito to the susceptible human 
population, and ^ the lifespan of an adult mosquito. 

If TZq < 1, on average, an infected individual produces less than one new infected individual 
over the course of its infectious period, and the disease declines. Conversely, if TZq > 1, then each 
individual infects more than one person, and the disease invades the population. 

Theorem 5. If A4 > and TZq > 1, then System ([T])-([2| admits the endemic equilibrium E^ — 
{S*h,rh,R*h,Al^,S;^,r*^) given by Eq. 0. 

Proof. The only solution of Eq. ([3| with Ih > or > 0, the only endemic equilibrium, is E^. 
That occurs, in agreement with Theorem [2j in the case > and x > The condition x > ^ 
is equivalent, by Theorem |4] to TZq > 1. □ 



Using the methods developed in Li and Muldowney ( 1996 1 and van den Driessche and Wat 



mough (2002 1, if 7?.o < 1, then the DFE is globally asymptotically stable in ft, and the vector-borne 
disease always dies out; if TZq > 1, then the unique endemic equilibrium is globally asymptotically 
stable in ft, so that the disease, if initially present, persists at the unique endemic equilibrium 
level. 



5 Numerical implementation 

In the epidemic of Cape Verde, infections were rising at a rate of one thousand people a day. 
The entire population was asked to participate in the campaign of cleaning and eradication of the 
mosquito, with the help of the police and the army. Data for humans are available at .www . ine . cv 
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(a) Human population in the compartments Susceptible (Sh), Infected (Ih), and Recovered {Rh)- 



gXlo Number of mosquitoes 
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Time (days) 



(b) Mosquito population in the compartments Aquatic phase (Am), Susceptible (Sm), and Infected (Im)- 

Figure 2: Variation of human and mosquito populations when no control is applied, i.e., for System 
with CA = 0, Cm = 0, and a = 1. 



(2011 ), but knowledge of mosquitoes in Africa is poor. For Aedes aegypti, Esteva and Yang ( 2005[ ) 



and Coelho et al. (2008) have collected observations from Brazil. The simulations were carried 
out with Nh = 480,000, B 0.8, /3„,i = 0.375, Phm = 0.375, i^h = 1/(71 x 365), rjt = 1/3, 
Mm = 1/10, /im = 1/10, if = 6, ila ~ 1/4, rjA = 0.08, to = 3, fc = 3. The initial conditions for the 
problem are: Sho = Nh - 10, ho = 10, Rm = 0, A„m = kNh, S,nO = mNh, Imo = 0. With these 
values, Eq. ^ gives > 0. The time interval is one year and tf — 365 days. We performed all 
simulations and graphics with Matlab. To solve System ([T])-([2]), we used the ode45 routine. This 
function implements a Runge-Kutta method with a variable time step for efficient computation. 

Figures [2a| and |2b| show the human and mosquito populations in the absence of any control. 
The human infection reached a peak between the 30th and the 50th day. The infection of the 
mosquitoes is delayed. The total number of infected humans obtained from System ([I])-([2]) is 
higher than observations in Cape Verde. The difference can be due to our inability to quantify 
individual prophylactic efforts. 

Figures [3a| and [3b| on adulticide control. Figures [4a| and [4b] on larvicide control, and Figures [5a| 
and |5b| on mechanical control show that a small quantity of each control is efficient to drop 
infection. Figures [3a| and [3b| show that the human population is already well protected by covering 
only 25% of the country with insecticide for adult mosquitoes. However, we consider that Aedes 



7 



xio" Infected Human 




Time (days) 

(a) Variation of infected individuals (//j) with adulticide. 




° ">° Time(days)^=° 

(b) Variation of infected mosquitoes (/m) with adulticide. 

Figure 3: Variation of infected human and infected mosquito populations with different levels of 
adulticide: Cm = 0, Cm = 0.25, Cm = 0.50, Cm = 0.75, and Cm = 1- The other controls are not 
applied: = 0, and a = 1. 
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10 20 30 40 50 60 70 80 90 100 

Time (days) 

(b) Variation of infected mosquitoes {Im.) with larvicide. 

Figure 4: Variation of infected human and infected mosquito populations with different levels of 
larvicide: ca = 0, ca = 0.25, ca = 0.50, ca = 0.75, and ca = 1- The other controls are not 
applied: Cm = 0, and a = 1. 
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„ X 1 0* Infected Human 
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(a) Variation of infected individuals (7^) with mechanical control. 
X 1 0* Infected Mosquito 




Time (days) 

(b) Variation of infected mosquitoes {Im) with mechanical control. 

Figure 5: Variation of infected human and infected mosquito populations with different levels of 
mechanical control: a = 0.01, a = 0.25, a = 0.5, a = 0.75, and a = 1. The other controls are not 
applied: = 0, and ca = 0- 
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50 100 150 

Time (days) 

(a) Variation of infected individuals (7^) with a combined use of the three controls. 




50 100 150 

Time (days) 

(b) Variation of infected mosquitoes (/m) with a combined use of the three controls. 

Figure 6: Variation of infected human and infected mosquito populations when using larvicide, 
adulticide, and mechanical control simultaneously: ca = Cm = 1 — a = 0, 0.01, 0.05, 0.1, 0.15. 



aegypti does not become resistant to insecticide and that there is no shortage of insecticide. 
Figures 4a 4b 5a| and [5b] show the controls applied in the aquatic phase of the mosquito. The 
controls were studied separately, but each one is closely related to the other. None of these controls 
is sufficient to drop the total number of infected humans to zero, but the removal of breeding sites 
and the use of larvicide contributes to prophylaxis. 

Figures |6a| and |6b| show simulations using the three controls simultaneously. They show that 
10% of each control, applied continuously, is enough to contain the infection near zero case. 

The eradication of Aedes aegypti may not be feasible and, from the environment point of view, 
not desirable. The aim is to reduce the mosquito density and increase the immunity on the humans. 
Population herd immunity can be reached by increasing the total number of resistant persons to 
the disease, which implies that these persons have been infected, or through vaccination. 

No commercially available clinical cure or vaccine is currently available for dengue, but efforts 



are to develop one (Blaney et al. 2007 Hombach et al. 20071. Effective vaccines have been pro 



duced against other flavivirus diseases such as yellow fever, Japanese encephalitis, and tick-borne 
encephalitis, so there is good hope for a vaccine for dengue. We now simulate the epidemiological 
model with vaccination. 



11 



Figure 7: Epidemiological SVIR (Susceptible, Vaccinated, Infected, Recovered) model for humans 
with Sh, Vh (human vaccinated population), Ih, and Rh as state variables, p (proportion of the 
vaccinated new born) as the control, and Nh, B, Pmh, fJ-h, Vh, w (proportional rate at which vacci- 
nation loses effect), 5* (fraction of the vaccinated susceptible), and a (infection rate of vaccinated 
people) as given parameters. 



6 Model with vaccination 



While direct individual protection is the major focus of mass vaccination program, population 
effects also contribute to individual protection through herd immunity, providing protection for 
unprotected individuals (Farrington 20031). The more vaccinated people, the less likely a suscep- 



tible person will come into contact with the infection. With the introduction of a vaccine, the SIR 
model related to the human population is augmented into the SVIR model of Figure [Tj where Vh 
represents the total number of vaccinated people. 

Vaccination is continuous with a constant proportion p of vaccinated new born. A fraction ^ 
of the susceptible is vaccinated. The vaccination reduces but does not eliminate susceptibility to 
infection. For this reason we consider the infection rate a of vaccinated people: when cr = the 
vaccine is perfect and when a = 1 the vaccine has no effect at all. The vaccination loses efficacy 
at a rate w. The differential system for the host population is: 

'S'hit) = il-p)^ihNh + wVh{t)- (^Bp,^h'^ + '^ + I^h) Shit) 
V{^{t) = p^XhNh + ^Sh{t) -(w + aB/3mh^ + f^h) 14(0 
/;(<) = B/3,r.h^{Shit) + aVhit)) ~ (j]h + \ih)lh{t) 
R'^{t)^^^hhit)-^^hRh{t). 

Figure [8] shows simulations for decreasing vaccine efficacy. Larvicide, insecticide, and mechan- 
ical control were kept null, and the parameters a and w were changed. For 80% of the population 
vaccinated, the efhcacy of the vaccine reduces the disease spread. 

Figure |9] presents the proportion of population vaccinated. It shows that dengue fever prophy- 
laxis articulates health and sustainable development. 



7 Conclusion 

Our simulations, based on our compartmental epidemiological model and formalizing clean-up 
campaigns to remove the vector breeding sites and the application of insecticides (larvicide and 
adulticide), have shown that even with a low continuous control over time, the results are encour- 
aging. The model with an imperfect vaccine has shown that the total number of infected persons 
can decrease quickly. 
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Figure 8: Total number of infected humans for decreasing vaccine efficacy {a = 1 — w = 0, 0.25, 
0.5, 0.75, 1), p = 0.80, and * = 0.80. 



xio" Infected Human 




50 100 150 

Time (days) 



Figure 9: Total number of infected humans when using different values of population vaccinated 
= 0, 0.25, 0.5, 0.75, l),p = 0,w = 0.85, and a = 0.15. 
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